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\0 ■ ABSTRACT 



Context. The predicted orbital period histogram of a subdwarf B (sdB) population is bimodal with a peak at short ( < 10 days) and 
long ( > 250 days) periods. Observationally, however, there are many short-period sdB systems known, but only very few long-period 

■ sdB binaries are identified. As these predictions are based on poorly understood binary interaction processes, it is of prime importance 
to confront the predictions to well constrained observational data. We therefore initiated a monitoring program to find and characterize 
long-period sdB stars. 

Aims. In this contribution we aim to determine the absolute dimensions of the long-period binary system PG 1 104+243 consisting of 
^ . an sdB and a main-sequence (MS) component, and determine its evolution history. 

5^ ■ Methods. High-resolution spectroscopy time-series were obtained with HERMES at the Mercator telescope at La Palma, and ana- 

■ lyzed to determine the radial velocities of both the sdB and MS components. Photometry from the literature was used to construct the 
_ spectral energy distribution (SED) of the binary. Atmosphere models were used to fit this SED and determine the surface gravity and 

temperature of both components. The gravitational redshift provided an independent confirmation of the surface gravity of the sdB 
component. 

Results. An orbital period of 753 ± 3 d and a mass ratio of ^ = 0.637 ± 0.015 were found for PG 1 104+243 from the radial velocity 
^ I curves. The sdB component has an effective temperature of T^fi = 33500 ± 1200 K and a surface gravity of log^ = 5.84 ± 0.08 dex, 

0^ ■ while the cool companion is found to be a G-type star with T^ff = 5930 ± 160 K and log g = 4.29 ± 0.05 dex. When a canonical mass 

■ of MsdB = 0.47 Mq is assumed, the MS component has a mass of Mms = 0.74 ± 0.07 Mg, and its temperature corresponds to what is 
0^ ' expected for a terminal age main-sequence star with sub-solar metalicity. 

Conclusions. PG 1 104+243 is the first long-period sdB binary in which accurate and consistent physical parameters of both compo- 
nents could be determined, and the first sdB binary in which the gravitational redshift is measured. Furthermore, PG 1104+243 is the 
first sdB+MS system that shows consistent evidence for being formed through stable Roche-lobe overflow. An analysis of a larger 
sample of long-period sdB binaries will allow for the refinement of several essential parameters in the current formation channels. 
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, 1 . Introduction leal phase, but at a later stage, when the companion penetrates 
^ ■ the initial radiative layer of the giant. These channels produce 
■ ■ ■ Subdwarf B (sdB) stars are considered to be core helium burning short-period sdB binaries with = 0. 1 - 10 d. Furthermore, 
stars with a very thin (Mh < 0.02 Mq) hydrogen envelope, and a there are two stable Roche-lobe overflow (RLOF) channels. In 
mass close to the core helium flash mass of 0.47 Mq. Their exis- the first, distinct mass transfer starts near the tip of the first gi- 
tence provides stellar evolutionai-y theory with a challenge: how ant branch, and in the second, mass transfer occurs when the 
can a red giant lose its entire hydrogen e nvelope before starting binary is passing through the Hertzsprung gap. Subdwaif B bi- 
core helium fusion?'Mengel et alJ was the first to propose naiies formed in this way are predicted to have longer periods 
the hypothesis of binary evolution to form sdB stars. Currendy ranging from 10 to over 500 days. All of the above scenai'ios re- 
there are five evolutionary channels that are thought to produce suit in a very small mass range for the sdB component, peaking 
sdB stars, all of which are binary evoluti on cha nnels where inter- at M^ab = 0.47 ± 0.05 M^. The last channel that can produce an 
action physics play a major role (He be3l2009l) . The first two are sdB star is the double white dwarf merger channel, where a pair 
common envelope (CE) ejection channels. In the first of these, a of close helium-core white dwaifs loses orbital energy through 
CE is formed by unstable mass transfer and expelled due to en- gravitational waves, resulting in a merger. This channelforms a 
ergy transfer from the shrinking orbit to the envelope. In the sec- single subdwarf star with a higher mass, up to 0.65 Mq (|Webbink| 
ond CE channel, the envelope is not ejected in the initial dynam- ISM)- 
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The binary interaction processes that form these sdB bi- 
naries are very complex, and many aspects are still not com- 
pletely understood. Formally, they are described using several 
parameters such as efficiency of envelope ejection, accretion ef- 
ficiency, physical interaction during common envelope phase, 
etc. Most of these are currently poorly constrained by obser- 
vations. The most common tool to constrain evolution paths is 
by performing a binary population synthesis (BPS) and com- 
paring the resulting population to an observed popula t ion fo r 
different values of t hese param eters. 'Ha n et all (l2002i l2003h . 
iNelemansI i20 1 Ot) and lClausen et al. (2012) have performed such 
BPS studies. The main results are that subdwarf progenitors with 
a low-mass companion will have unstable mass transfer on the 
RGB, during which their orbit will shrink, resulting in short- 
period sdB binaries. Subdwarf progenitors with a more massive 
companion will have expanded orbits, resu lting in long-period 
sdB binaries. Several observational studies (Maxted et al.ll200lt 



Table 1. The observing date (mid-time of exposure), exposure 
time and the initials of the observer of the PG 1 104-1-243 spectra 
observed with HERMES. 



[Morales-Rueda et all 12003; Copperwheat et al. 20 lib have fo- 
cused on short-period binaries, and found that about half of all 
sdB stars reside in short-period (Poib < 10 d) binaries. Longe r 
period sdB binaries have been observed (e.g. lGreen et al ] l200ll) . 
but no definite orbits ha ve been established. 

IClausen et al.l (12012 ) concluded that, although the currently 
known population of sdB binaries with a white dwarf (WD) or M 
dwarf companion can constrain some parameters used in the bi- 
nary evolution codes, many parameter remain open. Furthermore 
they concluded that sdB binaries with a main-sequence (MS) 
component will be able to provide a better understanding of the 
limiting mass ratio for dynamically stable mass transfer during 
the RLOF phase, the way that mass is lost to space, and the trans- 
fer of orbital energy to the envelope during a common envelope 
phase. The period distribution of sdB + MS binaries can support 
the existence of stable mass transfer on the RGB. If long -period 
sdB + MS binaries are found, they are thought to be formed 
during a stable mass transfer phase in the RLOF evolutionary 
channel. If no long-period sdB + MS binaries can be found, this 
would indicate that mass transfer on the RGB is unstable. 

The goal of this paper is to describe the methods we used to 
determine radial velocities of both the main-sequence and sub- 
dwarf B component in sdB + MS systems (Sect.|2]i, and the 
use of spectral energy distributions (SEDs) to determine the 
spectral type of both components (Sect.|3]l. Furthermore, if the 
mass ratio is known from spectroscopy, the SEDs from high- 
precision broad-band photometry alone can be used to deter- 
mine their surface gravities with an accuracy up to Alog^ - 
0.05 dex (cgs) and the temperatures with an accuracy as high as 
5%. The resulting surface gravities are independently confirmed 
by the surface gravities derived from the gravitational redshift 
(Sect.lHi. After which all obtained parameters are summarized 
(Sect.jSj. The equivalent widths of iron lines in the spectra can 
be used to independently confirm the atmospheric parameters 
of the cool companion (Sect.|6]l. These methods were applied 
to PG 1 104+243. The cool companion of PG 1 104+243 is pre- 
sumed to be on the main-sequence, and we will refer to it as 
the MS component in this paper PG 1 104+243 is part of a long- 
term spectroscopic observing program, and preliminary results 
of this and seven more systems in this pr ogram were presented 
bv l0stensen & Van Winckel (1201 iL 1201 2h . 

2. Spectroscopy 

The high-resolution spectroscopy of PG 1104+243 was ob- 
tained with the HERMES spectrograph (High Efficiency and 
Resolution Mercator Echelle Spectrograph, R = 85000, 55 
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orders, 3770-9000 A, iRaskin et al.l 1201 ih attached to the 
1.2-m Mercator telescope at the Roque de los Muchachos 
Observatory, La Palma. HERMES is isolated in an over pressur- 
ized temperature-controlled enclosure for optimal wavelength 
stability, and connected to the Mercator telescope with an op- 
tical fiber. In total 38 spectra of PG 1 104+243 were obtained be- 
tween January 2010 and February 2012. The date and exposure 
time of each spectrum is shown in Table[T] HERMES was used 
in high-resolution mode, and Th-Ar-Ne exposures were made at 
the beginning and end of the night, with the exposure taken clos- 
est in time used to calibrate the wavelength scale. The exposure 
time of the observations was adapted to reach a signal-to-noise 
ratio (S/N) of 25 in the V-band. The HERMES pipeline v3.0 was 
used for the basic reduction of the spectra, and includes bary cen- 
tric correction. Part of a sample spectrum taken with HERMES 
is shown in Fig.[T] 

Furthermore, five intermediate-resolution spectra (R - 4100) 
were obtained with the Blue Spectrograph at the Multi Mirror 
Telescope (MMT) in 1996-1997. The Blue Spectrograph was 
used with the 832/mm grating in 2nd order, covering the wave- 
length region 4000-4950 A. The data were reduced using stan- 
dard IRAF tasks. An overview of the spectra is given in Table[3] 

To check the wavelength stability of HERMES, 38 differ- 
ent radial velocity standard stars of the lAU were observed over 
a period of 1481 days, one standard star a night, coinciding 
with the observing period of PGl 104+243. These spectra are 
cross-correlated with a line mask corresponding to the spec- 
tral type of each star to determine the radial velocity. For this 
cross-correlation (CC), the hermesVR method of the HERMES 
pipeline, which is designed to determine the radial velocity of 
single stars, is used. This method handles each line separately, 
starting from a given line-mask. The method uses the extracted 
spectra, after the cosmic clipping was performed but prior to 
normalization. It performs a cross-correlation for each line in 
the mask, and sums the CC functions in each order The ra- 
dial velocity is derived by fitting a Gaussian function to the CC 
function (Raskin et al. 2011). To derive the final radial velocity 
only orders 55— >74 (5966 - 8920 A) are used, as they give the 
best compromise between maximum S/N for G-K type stars and. 
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Fig. 1. A sample normalized spectrum of PG 1104+243 (black), showing the 5875 A He I line of the sdB component, and several 
lines as e.g. the Na doublet at 5890-5896 A of the MS component. Furthermore, a synthetic spectrum of the MS component (green) 
and sdB component (red) are shown. 



after maskin g the t elluric bands, absence of telluric influence 
dRaskin et al.il201 fl) . The resulting standard deviation of these 
radial velocity measurements is 80 m/s with a non significant 
shift to the lAU radial velocity standard scale. 

2.1. Radial Velocities 

The determination of the radial velocities of PG 1 104H-243 from 
the HERMES spectra was performed in several steps. First, there 
were several spectra taken in a short period of time. As the 
or bital period of PG 1104-1-243 wa s determined at 752 + 14 d 
by I0stensen & Van Winckel ( 1201 2'). spectra that were taken on 
the same night, or with only a couple of nights in between are 
summed to increase the S/N. We experimented with several dif- 
ferent intervals, and determined that when spectra taken within 
a five-day interval (corresponding to 0.7% of the orbital period, 
or a maximum radial velocity shift of 0.05 km/s) were merged, 
there was no significant smearing or broadening of the spectral 
lines. After this merging, 25 spectra with a S/N ranging from 25 
to 50 remained. 

The MS component has many spectral lines in every Echelle 
order, most of which are not disturbed by the He lines of the sdB 
component. Because of this, the CC-method of the HERMES 
pipeline (hermesVR) can be used to determine the radial veloc- 
ities of the MS component without any loss of accuracy. For 
the MS component of PG 1 104-1-243, a G2-type mask was used. 
The final errors on the radial velocities take into account the for- 
mal errors on the Gaussian fit to the normalized cross-correlation 
function, the error due to the stability of the wavelength calibra- 
tion and the error arising due to the used mask. The resulting 
velocities and their errors are shown in Table|2] 

The determination of the radial velocities of the sdB compo- 
nent is more complicated. Except for the Balmer lines, there are 
only a few broad He lines visible in the spectrum. After compar- 
ison of a synthetic spectrum of a G2-type star with the spectra, 
it was found that only the 5875.61 A He I blend is not contam- 
inated by lines of the MS component. For this single blend, we 
had to deploy another method to obtain an accurate radial veloc- 
ity. As only one line is used, anomalies in the spectra can cause 
significant errors in the radial velocity determination. To prevent 
this the region around the He I line is cleaned from any remain- 



ing cosmic rays by hand. The normalization process is designed 
to run automatically, making it possible to quickly process mul- 
tiple lines is many spectra, and consists of two steps. In the first, 
a spline function is fitted through a 100 A region centered at 
the He I line to remove the response curve of HERMES, and the 
main features of the continuum. In the second step, a small 10 A 
region around the He I line is further normalized by fitting and 
re-fitting a low degree polynomial through the spectrum. After 
each iteration, the flux points lower than the polynomial are dis- 
carded. This is repeated until the polynomial fits the middle of 
the continuum. An example of a spectrum normalized using this 
method is shown in Fig.[T] 

These cleaned and normalized spectra can then be cross- 
correlated to a synthetic sdB spectrum. For this purpose a high- 
resolution LTE_s£ectrim ofTgff = 35000 K and \ogg = 5.50 from 
the grids o f iHeber et all (l2000l) was used. The resolution of this 
synthetic spectrum is matched to the resolution of the HERMES 
spectra. Synthetic spectra with different temperatures and sur- 
face gravities were tested, to check the temperature and \ogg 
dependence of the radial velocities. Changing either of those pa- 
rameters causes a systemic shift of maximum 0.2 km/s in the 
resulting radial velocities. The cross-correlation is performed 
in wavelength space by starting at the expected radial velocity, 
calculated from the radial velocities of the MS component. To 
determine the radial velocity, a Gaussian is fitted to the cross- 
correlation function. The error on the radial velocities is obtained 
from a Monte-Carlo simulation. This is done by adding Gaussian 
noise to the spectrum and repeating the cross-correlation with 
the synthetic spectrum. The level of the Gaussian noise is deter- 
mined from the noise level in the continuum near the He I line. 
The final error is determined by the standard deviation from the 
radial velocity results of 1000 simulations, the wavelength sta- 
bility of HERMES and the dependence on the template. The re- 
sulting radial velocities and their errors are shown in Table|2] 

To derive the radial velocities for the MS component from 
the MMT spectra, the spectra were cross-correlated against two 
high S/N GOV specti-a (HD 13974 and HD39587) using the IRAF 
task FXCOR, where the cross-correlation was done only in the 
wavelength ranges containing absorption lines from the MS 
companion, ie avoiding the Balmer lines. He A 4026, He A 4387, 
He^ 4471, He^ 4686, He A 4712, and He^ 4921. The veloc- 
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Table 2. The radial velocities of both components of 
PG 1104+243. 
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Table 3. The observing date (mid-time of exposure), ex- 
posure time and radial velocities of the MS component of 
PG 1104-H243, determined from the spectra observed with the 
Blue Spectrograph at the MMT. 
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ities derived from the two GOV stars were averaged to get the 
velocity of the MS companion for each of the five PG 1 104-1-243 
spectra. The proper scaling factor for the velocity eiTors output 
by FXCOR was determined using the standard deviation of the ve- 
locities obtained from each template relative to the average ve- 
locity for each epoch. The resulting radial velocities and errors 
and are given in Table[3] 

2.2. Orbital Parameters 

The orbital parameters were calculated by fitting a Keplerian or- 
bit through the radial-velocity measurements, while adjusting 
the period (P), time of periastron (To), eccentricity (e), angle 
of periastron (cj), amplitude (K) and systemic velocity (y), us- 
ing th e orbital period determined by I0stensen & Van Winckell 
(12012 *) as a first guess for the period. The MS and sdB compo- 
nent were treated separately in this procedure. The radial veloci- 
ties of both components were weighted according to their errors 
as w = 1 /cr. Fitting the orbit resulted in a very low eccentricity 



Table 4. Spectroscopic orbital solution for both the main- 
sequence (MS) and subdwarf B (sdB) component of 
PG1104-H243. 
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Notes, a denotes the semi-major-axis of the orbit. The quoted errors are 
the standard deviation from the results of 5000 iterations in a Monte 
Carlo simulation. 



(e < 0.002) . To check if the orbit of the sdB+MS binary is circu- 
larized, the ILucv & Sweenev. (.1971,1 test for circularized orbits 
was applied: 



(1) 



where ecc indicates the residuals of an eccentric fit, and circ the 
residuals of a circular fit, n the total number of observations, 
and m - 6, the number of free parameters in an eccentric fit. If 
P < 0.05, the orbit is significantly eccentric, if f > 0.05, the 
orbit is effectively circular In the case of PG 1 104-1-243 we find 
P - 0.15, indicating a circular orbit. This is expected as the orbit 
is supposed to be circulariz ed by tidal interactions during the 
RGB evolution (IZahnlll977h . In the further determination of the 
orbital parameters, a circular orbit is assumed. 

To derive the final orbital parameters, both components were 
first treated separately to measure the difference in systemic ve- 
locity (see Sect.|4]i. Then that difference was subtracted from the 
radial velocities of the sdB component, and both components are 
solved together to obtain their spectroscopic parameters. The un- 
certainties on the orbital parameters were determined by using 
Monte Carlo-simulations. The radial velocities were perturbed 
based on their errors, and the errors on the parameters were de- 
termined by their standard deviation after 5000 iterations. The 
residuals of the orbital fit to the MS component are larger than 
the errors on the radial velocities. However, this variability is 
most likely intrinsic, caused by stellar spots or low amplitude 
pulsations, and does not affect the derived orbital parameters. 
The spectroscopic parameters of PG 1104-1-243 are shown in 
TablelH The radial-velocity curves and the best fit solution are 
plotted in Fig|2l 



3. Spectral energy distribution 

To derive the spectral type of the main-sequence and subdwarf 
component, we fitted the photometric spectral energy distribu- 
tion (SED) of PG 1 104+243 with model SEDs. With this method 
we can determine the effective temperature and surface gravity 
of both components with good accuracy. 
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Fig. 2. Top: spectroscopic orbital solution for PG 1104+243 
(solid line: MS, dashed line: sdB), and the observed radial ve- 
locities (blue filled circles: MS HERMES spectra, blue open cir- 
cles: sdB HERMES spectra, green triangles: MS MMT spectra). 
The measured system velocity of both components is shown by 
a dotted line. Middle: residuals of the MS component. Bottom: 
residuals of the sdB component. 



3.1. Photometry 

Photometr y of PG 1 104-H2 43 was collected using the subdwarf 
databas43 (0stens enl2OO6h . which contains a compilation of data 
on hot subdwarf stars collected from the literature. In total, fif- 
teen photometric observations were found in four different sys- 
tems: Johnson B and V, Cousins R and /, Stromgren uvby, and 
2MASS J, H and Kj. Four observations have uncertainties larger 
or equal than 0. 1 mag and are discarded. The observations used 
in the SED fitting process are shown in Table|5] Accurate pho- 
tometry at both short and long wavelengths are used to establish 
the contribution of both the hot sdB component and the cooler 
MS component. 



3.2. SED fitting 

In the SED of PG 1 104-H243 in Fig.H we see both the steep rise 
in flux towards the shorter wavelengths of the sdB component 
and the bulk in flux in the red part of the SED of the MS compo- 
nent. To fit a synthet ic SED to the observed photometry, Kurucz 
atmosphere models (lKuruczlll993h are used for the MS compo- 
nent, and TM AP (Tubingen NLTE Model-Atmosphere Package, 
IWerne r et al.'^2003) atmosphere models for the sdB component. 
The Kurucz models used in the SED fit have a temperature range 
from 4000 to 9000 K, and a surface gravity range of log ^=3.5 
dex (cgs) to 5.0 dex (cgs). The TMAP models we have used 
cover a temperature range from 20000 K to 50000 K, and log g 
from 4.5 dex (cgs) to 6.5 dex (cgs). They are calculated using 
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Table 5. Photometry of PG 1 104+243 collected from the litera- 
ture. 



Band 


Wavelength 


Width 


Magnitude 


Error 




A 


A 


map 


map 


Johnson B" 


4450 


940 


11.368 


0.010 


Johnson V" 


5500 


880 


11.295 


0.013 


Cousins R'^ 


OJVU 


1 Q on 


11 1 A 1 

1 i. ioi 


n n 1 n 
U.UlU 


Cousins /" 


7880 


1490 


11.001 


0.010 


Stromgren m* 


3460 


300 


11.513 


0.045 


Stromgren v* 


4100 


190 


11.485 


0.045 


Stromgren fc* 


4670 


180 


11.334 


0.045 


Stromgren ;y* 


5480 


230 


11.256 


0.006 


2MASS r 


12410 


1500 


10.768 


0.026 


2MASS 


16500 


2400 


10.520 


0.027 


2MASS K,' 


21910 


2500 


10.510 


0.023 



R eferences. lAUard et alj ( 119941) 

''^> ISkrutskieet"an (120061) 



IWesemael et"aLl ( Il992h 



an atmospheric mixture of 97 % hydrogen and 3 % helium in 
mass due to the expected He-depletion in the atmosphere of 
sdB stars. The original TMAP atmosphere models cover a wave- 
length range of 2500-15000 A. To include the 2MASS photome- 
try as well, they are extended with a black body of corresponding 
temperature to 24000 A. 

The S EDs are fitted f ollow ing the grid-based method de- 
scribed in Degroote et alJ (1201 ih . but we extended it to include 
constraints from binarity. In the binary scenario, there are eight 
free parameters to consider; the effective temperatures (Teff MS 
and Teff^sdE), surface gravities (gus and gsdfi) and radii {Rms and 
^sdfi) of both components. The interstellar reddening E(B - V) is 
naturally presumed equal for both components, and is incorpo- 
rated using the reddening law of Fitzpatrick (2004). The models 
are first corrected for interstellar reddening, and then integrated 
over the photometric passbands. The distance (d) to the system, 
acts as a global scaling factor. 

The radii of both components are necessary as scale factors 
for the individual fluxes when combining the atmosphere models 
of both components to a single SED. Including the distance to 
the system, the flux of the combined SED is given by: 



Fiot(A) - (^MS^Ms('t) + ^sdB-^sdB('t)) , 



(2) 



where Fms is the flux of the MS model, and FsdB is the flux of 
the sdB model, in a particular passband. However, if the masses 
of both components are known, the radii can be derived and re- 
moved as free parameters. Assuming that the mass of the sdB 
component is /V/sdB - 0.47 Mq as predicted by stellar evolu- 
tion models (see Sect. [1]), the mass of the MS component can 
be calculated using the mass ratio derived from the radial ve- 
locity curves (see Sect. 12.2b : /Wms = 1/'? A^sdB- For each model 
in the grid, the radii of both components are derived from thek 
surface gravities and masses, according to: 



GMi 



, / = sdB, MS, 



(3) 



in which G is the gravitational constant. This relation can be used 
to eliminate the radius dependence in Eq.|2| 



gM5 \ gsdB 



(4) 
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The distance to the system is computed by shifting the combined 
synthetic model flux {F,„,) to the photometric observations. As 
the effect of the mass of the main-sequence component (Mms) in 
this equation can be adjusted by the distance, the resulting flux 
is only dependent on the mass ratio of the components, and not 
the presumed mass of the sdB component. By using this mass 
ratio as a limiting factor on the radii of both components, the 
number of free parameters in the SED fitting process is reduced 
from eight to six. 

To select the best model, the^^ value of each fit is calculated 
using the sum of the squared errors weighted by the uncertainties 
(^,) on the observations: 



X 



{Oi - Cif 



(5) 



with O, the observed photometry and C, the calculated model 
photometry. The expectation value of this distribution is ^ = 
Nobs - Nfree, with Nobs the number of observations and Nfree the 
number of free parameters in the fit. In our case, k - 11-6 = 5. 
Based on this statistics the error bars on the photometry can 
be checked. If the obtained is much higher than the expecta- 
tion value, the errors on the photometry are underestimated. If 
on the other hand, the obtained is much lower than the expec- 
tation value, the errors on the photometry are overestimated or 
we are over fitting the photometry. 

If the ;if^-statistics are valid, the uncertainties on the fit pa- 
rameters can be estimated using the cumulative density function 
(CDF) of the statistics. To calculate the CDF, the obtained 
values have to be rescaled. The of the best fitting model is 
shifted to the expectation value of the ;t'^-distribution, and all 
other values in the grid are scaled relative to the h^'sXx^. The 
probability of a model to obtain a certain value is given by 
the CDF, which can be computed using; 



(6) 



Where P is the regularized F-function. Based on this CDF the 
uncertainties on the parameters can be calculated based on the 
distribution of the probabilities. 

3.3. Resu/fs 

To fit the SED of PG 1 104-1-243, first a grid of composite binary 
spectra is calculated for ~2 000 000 points randomly distributed 
in the Tetf , log g and E(B-V) intervals of both components. These 
intervals are then adjusted based on the 95% probability intervals 
of the first fit, and the fitting process is repeated. To determine 
final confidence intervals, separate grids of ~ 1 000 000 points 
which only vary two parameters, while keeping the other pa- 
rameters fixed at their best fit values, are used. The uncertainties 
on the parameters are determined based on the 95% probabil- 
ity intervals. The best fit has a reduced of 6.3, and with five 
degrees of freedom, the values in the grid were rescaled by 
a factor 1 .26 which is equivalent to slightly increasing the pho- 
tometric uncertainties. The SED fit resulted in an effective tem- 
perature of 5930 + 160 K and 33500 + 1200 K for the MS and 
sdB components respectively, while the surface gravity was de- 
termined at respectively 4.29 + 0.05 dex and 5.81 + 0.05 dex for 
the MS and sdB component. The reddening of the system was 
found to be 0.001^[!I!,'J mag . Thi s redde ning can be compared 
to the dust map 

of SchlegeleFiD dHH, which gives an upper 
bound of E(B-V) - 0.018 for the reddening in the direction of 




3 ■ 103 



1 ■ 10* 
Wavelength (A) 



Fig. 3. Top: The spectral energy distribution of PG 1104H-243. 
The measurements are given in blue, the integrated synthetic 
models are shown in black, where a horizontal error bar indi- 
cates the width of the passband. The best fitting model is plotted 
in red. Bottom: The (O-C) value for each synthetic flux point. 



PG 1104+243. The reddening of E(B-V) = 0.001+° °^^ mag in- 
dicates that PG 1104-1-243 is located in front of the interstellar 
molecular clouds. This is confirmed spectroscopically, as there 
are no sharp interstellar absorption lines (e.g. Ca n K-H, K i) 
visible in the spectra. These parameters indicate that the cool 
companion is a G-type star 

The results of the fit together with the probability intervals 
are shown in Table|6] The optimal SED fit is plotted in Fig. [3] 
while the grids with the probability for each grid point are plot- 
ted in Fig.|4] The probability distribution for the surface grav- 
ity and effective temperature of the MS component follows a 
Gaussian pattern, while the probability distribution for the same 
parameters of the sdB component is elongated. This latter ef- 
fect is related to a moderate correlation between logg and Teff. 
Thus, the effect on the models caused by an increase in Teff , can 
be diminished by increasing logg as well. This correlation is 
stronger towards shorter wavelengths, so that the main-sequence 
component, which is mainly visible in the red part of the SED, 
is less affected. A similar effect is visible in the relation between 
E{B - y) and Teff. 

As a check of the results, a SED fit of PG 1 104+243 without 
any assumptions on the mass of both components is performed 
as well. In this case, the radii of the components are randomly 
varied between 0.5 and 2.0 for the MS component and be- 
tween 0.05 /?0 and 0.5 /?o for the sdB component. The results of 
this fit are given in the lower half of Table|6] The resulting best 
fit parameters are very close to the results from the mass ratio- 
constrained fit, but apart from the uncertainty on the MS effec- 
tive temperature, the uncertainties on the parameters are much 
larger. For both the MS and sdB surface gravity, the 95% prob- 
ability interval are larger than the range of the models. The ef- 
fective temperature and surface gravity of both components of 
PG 1 104+243 correspond with the ionization balance of the iron 
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Fig. 4. The confidence intervals of the SED fit of PG 1 104+243. The upper row show distributions for the main-sequence compo- 
nent, the lower row for the subdwarf component. The best fit results are indicated with a red cross. The different colors show the 
cumulative density probability connected to the;^'^ statistics given by equation|6] with the;i'^ values of the grids rescaled by a factor 
of 1.26. 



lines seen in the spectrum of MS component, and the B aimer 
lines of the sdB component. 



The reason for the high accuracy of the mass ratio- 
constrained fit compared to a fit in which the radii are uncon- 
strained is shown in Fig.|5] Here five model SEDs are shown, 
calculated using the best fit parameters from Table|6] but with 
varying surface gravity, and accordingly adjusted radii, for the 
sdB component. The best fit model with log^ = 5.81 dex (cgs) 
is shown in blue, while models with a log g at the edge of the 
confidence interval {\ogg - 5.77 and 5.85 dex) are plotted in 
black and green. When varying log g, the radius of the sdB com- 
ponent changes, while the radius of the MS component remains 
constant. It is this radius change that causes the large deviation 
in the models, and allows for an accurate determination of the 
surface gravity. When only the gravity is changed, and the radii 
of both components is kept constant, the difference in absolute 
magnitude in the Stromgren-u band is of the order of 10 mag, 
while when the radius is adapted using the mass, the difference 
is on the order of 10"^ mag. Furthermore, because the absolute 
flux depends on the radius, this method has the advantage that all 
photometric observations are used to constrain the log g of each 
component, instead of only the observations in the blue or the 
red part. This analysis of PG 1 104-1-243 clearly shows the power 
of binary SED fitting when the mass ratio of the components is 
known. 



4. Gravitational Redshift 

In an sdB + MS binary the difference in surface gravity between 
the two components is substantial. The surface gravity of a star 
gives rise to a frequency shift in the emitted radiation, which is 
known as the gravitational redshift. General relativity shows that 
the gravitational redshift as a function of the mass and surface 
gravity of the star is given by (lEinsteinlll9I6l) : 

= ^ VCMi- (7) 

Where Zg is the gravitational redshift, c the speed of light, G 
the gravitational constant, M the mass, and g the surface gravity. 
This Zg will effectively change the apparent systemic velocity for 
the star. In a binary system the difference in surface gravity for 
both components will be visible as a difference in systemic ve- 
locity between both components. As the Zg is proportional to the 
square root of the surface gravity, this effect is only substantial 
when there is a large difference in log^ between both compo- 
nents, as is the case for compact subdwarfs and main-sequence 
stars. 

The measured difference in systemic velocity can be used to 
derive an estimate of the surface gravity of the sdB component. 
Using the mass ratio from the spectroscopic orbit, and the canon- 
ical mass of the sdB component, the mass of the MS component 
can be calculated. Combined with the surface gravity of the MS 
component from the SED fit, the Zg of the MS component can 
be calculated. The Zg of the sdB component can be obtained by 
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Table 6. The results of the SED fit, together with the 95% and 
80% probability intervals. In the upper part (fixed radius) the 
results when the radii of both components are determined from 
the log g and mass, are shown. In the lower part (free radius) the 
results with unconstrained radii are shown. For some parameters 
the 95% and 80% probability intervals could not be determined 
as they are larger than the range of the models. 



Parameter 



Best fit 



95% 



80% 



Fixed radius 
Main sequence component 
reff(K) 5931 5769-6095 

logg(dex) 4.29 4.26-4.32 

E(B-V) (mag) 0.001 0.000-0.012 



Subdwarf B component 
T,ff (K) 33520 
logg(dex) 5.81 
E(B-V) (mag) 0.001 



32400-34800 

5.77-5.85 
0.000-0.013 



Free radius 
Main sequence component 
Teff (K) 5970 5730-6190 

logg(dex) 4.38 / 

E(B-V)(mag) 0.001 0.000-0.100 



Subdwarf B component 
r,ff(K) 33810 
logg(dex) 5.90 
E(B-V) (mag) 0.001 



27500-45100 

4.70-/ 
0.000-0.100 



5821-6035 
4.27-4.31 
0.000-0.008 



32764-34290 

5.79-5.83 
0.000-0.008 



5840-6130 
3.55-/ 
0.000-0.090 



29000-40100 

5.25-/ 
0.000-0.090 




3600 4000 4400 

Wavelength (A) 



4800 



Fig. 5. SED models calculated for the best fit parameters of 
PG 1 104+243 as given in Table|6] but with varying surface grav- 
ity of the sdB component. From bottom to top, the models have 
a log ^ of respectively 6.00, 5.85, 5.81, 5.77 and 5.00 dex (cgs). 
The crosses show the integrated model photometry. The photo- 
metric observations with their error bars are given in red. 



combining the Zg of the MS component with the measured dif- 
ference in systemic velocity, and can then be converted to an 
estimated surface gravity of the sdB component. 

Using a canonical value of 0.47 Mq for the sdB component, 
and the mass ratio derived in Sect. l2.2l we find a mass of 0.74 
+ 0.07 Mq for the MS component. The SED fit resulted in a 
log^MS of 4.29 + 0.05 dex. Using equation |7l the gravitational 
redshift of the MS component is czg = 0.46 + 0.05 km s"'. The 
measured difference in systemic velocity is 1.97 + 0.18 km s"', 
resulting in a total gravitational redshift of cZg = 2.43 ± 0.20 km 
s"' for the sdB component. This is equivalent to a surface gravity 
of loggsdB - 5.90 ± 0.08 dex, which is very close to the results 
of the SED fit. 

Other phenomena that can affect the presumed redshift of the 
MS compone nt, as for example conv ectional blueshift (see, e.g., 
Ibravins 1982i:lTakeda & Uendl2012h . are not taken into account 
as the contribution of such effects is predicted to be around the 
one-sigma level. 



5. Absolute dimensions 

Combining the results from the SED fit with the orbital param- 
eters derived from the radial velocities, the absolute dimensions 
of PG 1 104+243 can be determined. With an assumed mass for 
the sdB component, the inclination of the system can be derived 
from the reduced mass determined in Sect. 12.21 Using this incli- 
nation, the semi-major axis of the system can be calculated. The 
surface gravity of the sdB component as determined by the SED 
fit and the gravitational redshift correspond within errors. For 
the final surface gravity, the average of both values is taken. The 
systemic velocity of PG 1104+243 is determined by correcting 
the systemic velocity of the MS component with its gravitational 
redshift, resulting in y = -15.17 + 0.07 km s"'. The radius of 
both components is derived from their mass and surface grav- 
ity. The absolute dimensions of PG 1104+243 are summarized 
in TableE] 

For an sdB star of the canonical mass and a surface gravity 
of log ^ = 5.85 dex a temperature around 33000 K is expected, 
consistent with the result of the SED fit and the strength of the 
He i/He ii lines. This high surface gravity and temperature is con- 
sistent with the binary population models of Brown (2008) that 
find that low metalicity produces sdB stars located at the high 
log g end of the extreme horizontal branch. The canonical sdB 
mass assumption together with the derived mass ratio implies 
that the MS component must have a sub-solar metalicity. At so- 
lar metalicity it should have a radius of 0.8 Rq and a temperature 
of ~4800K, which is clearly ruled out by the SED fit. If we drop 
the canonical sdB mass assumption, and force the MS star to 
have normal metalicity, both the MS and sdB star must be more 
massive in order to be consistent with the temperature of the MS 
component an d the spectros copic mass ratio. From the He-MS 
models of (Paczv ns5ll97lh . one finds that a core-He burning 
model with M > 0.6 Mq would have Teff > 40000K, which is 
ruled out by both SED fit and spe ctroscopy. Recent evolution- 
ary models bv lBerteUi et al.l (12008 *) for main sequence stars with 
various metalicities, are perfectly consistent with the parameters 
from the SED fit for an 0.82 Mq model, provided that the met- 
alicity is assumed to be substantially sub-solar, at around Z = 
0.005. Such tracks are shown in Fig.|6] where one can see that 
an 0.82 Mq star has an expected temperature of around 5900 K, 
and a surface gravity of log § = 4.30 near the end of the main se- 
quence, corresponding with an age of 15.0 + 2.5 Gyr for the MS 
component. This is similar to the Hubble time and given all un- 
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Table 7. Fundamental properties for both the main-sequence 
(MS) and subdwarf (sdB) component of PG 1 104+243. 



Systemic parameters 

P{d) 

To (HJD) 



7 (km s ') 
1 

a (Re) 

in 

E(B - V) 
rf(pc) 

Component parameters 
MS 



753 ±3 
2450386 + 8 

< 0.002 
-15.17 ±0.07 
0.64 ± 0.02 
322 + 12 
32 ±3 
0.001 ± 0.017 
305 + 10 



A-(km s-') 
M(Mo) 
logg (cgs) 
^(Ro) 
T^effCK) 
L(U) 
Vq (mag) 
My (mag) 



4.42 ± 
0.74 4 
4.29 i 
1.02 ± 
5930: 
1.15 ± 
12.06 ± 
4.64 4 



0.08 
0.07 
0.05 
0.06 
: 160 
0.13 
0.05 
0.05 



SdB 

6.9 ± 0.3 
0.47 ± 0.05" 
5.85 ± 0.08 
0.13 ± 0.02 
33500 ± 1200 
22.5 ± 3.5 
12.03 ± 0.08 
4.61 ± 0.10 



Notes. Assumed value based on evolutionary scenario. 

certainties in the models, indicates that the system is extremely 
old. 

The luminosity of both components can be calculated using 
L = Ana-R^T*, resulting in Lms = 1-15 ± 0.13 Lq and LsdB = 22.5 
+ 3.5 Lq. The apparent V magnitudes of both components are 
obtained directly from the SED fitting procedure and are given 
by: VcMS = 12.06 + 0.05 mag and Vo.sdB = 12.03 + 0.08 mag. 
The absolute magnitude can be obtained by integrating the best 
fit model SEDs over the Johnson V band, and scaling the result- 
ing flux to a distance of 10 pc, resulting in Mvms - 4.64 + 0.05 
mag and My^dB = 4.61 ± 0. 1 mag. The distance to the system can 
then be calculated from log d - (myMy + 5)/5, which places the 
system at a distance of d - 305 ± 10 pc. Obviously, the distance 
for both components comes out as exactly the same, as this is a 
fixed condition in the SED fitting process. 

The p roper motion of PG 1104-1-243 as measured by 



ine proper mc 
lH0g et all (120001) is: 



(^la,^ls) = (-63.4, -22.1) + (1.8, 1.9) mas yr"' (8) 

Using the method of iJohnson & SoderblomI (Il987h . these num- 
bers together with the measured value of y, can be used to com- 
pute the galactic space velocity vector 

(U, V, W) = (-61.8, -55.3, -52.5) + (2.8, 2.7, 1.3) km s"' (9) 

where U is defined as positive towards the ga lactic center. Using 
the va lues for the local standard of rest from lDehnen & BinnevI 
(119981) . we get 

(U, V, W)lsr = (-51.8, -50.1, -45.3) + (2.8, 2.8, 1.4) km s^HlO) 

Following the selection criteria of iReddv et alJ (l2006h . 
PG 1104+243 is bound to the galaxy, and belongs to the thick 
disk population. 

6. Atmospheric parameters from Fei-ii lines 

In this section the equivalent widths of iron lines are used to de- 
rive the atmospheric parameters and metalicity of the MS com- 
ponent. This method is independent of the sdB mass assumption 
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Fig. 6. Evolutionary tracks calculated by i Bertelli et alj (l2008h 
with a metalicity of Z = 0.005, for three different masses: 0.60 
Mq, 0.82 Mo and 1.00 Mq. The results obtained for the cool 
companion of PG 1 104+243 are shown in red. 



used in the SED fit, and is used to confirm the effective tempera- 
ture and surface gravity obtained with that method. Furthermore, 
the mass of the sdB component can be derived by using evolu- 
tionary tracks to fit the mass of the MS component, and the mass 
ratio from the radial velocity curves. 

Before the EWs can be measured, the continuum contri- 
bution of the sdB component must be subtracted from the 
HERMES spectra, which are then shifted to the zero velocity 
based on the MS radial velocities and suiiuned. The LTE abun- 
dance calculation routine MOOG (ISnedenll973h was used to de- 
termine the atmospheric parameters of the MS component using 
the EWs of 45 Fei and 10 Feii lines in the wavelength range 
4400-6850 A, excluding blended lines, and wavelength ranges 
contaminated by Balmer lines, Hei or Hen lines from the sdB 
component. The atmospheric parameter determination also pro- 
vides the metalicity of the MS component. The at omic data for 
the ir on lines was taken from the VALD line lists (iKupka et alJ 
1999'), and the used oscillator strengths were laboratory values. 
The EWs are measured via integration and abundances are com- 
puted by an iterative process in MOOG in which for a given 
abundance, the theoretical EWs of single lines are computed 
and matched to the observed EWs. The effective temperature 
is derived by the assumption that the abundance of individual 
Fe I lines is independent of lower excitation potential, only Fe 
I lines are used because they cover a large range in lower ex- 
citation potential, resulting in T^s - 6000 + 250 K. A surface 
gravity of logg = 4.5 ± 0.5 dex is derived by assuming ioniza- 
tion balance between the iron abundance of individual Fe i and 
Fe II lines. Assuming an independence between the iron abun- 
dance and the reduced EW, a microturbulence velocity of 2.0 
+ 0.5 km/s can be derived. The final metalicity is based on the 
average abundance of all used lines, and result in a metalicity 
of [Fe/H] = -0.58 ±0.11 dex, corresponding to Z = 0.005 ± 
0.001 using the conversion of Bertelli et al. ( 1994). The error on 
the metalicity is the line-to-line scatter for the Fe abundances 
for the derived atmospheric parameters. The used Fe I lines are 
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Table 8. Atmospheric parameters obtained based on the EW of 
iron lines, and derived properties. See Sect.|6]for further expla- 
nation. 





MS 


sdB 


M(Mo) 


0.86 + 0.15 


0.54 ±0.10 


logg (cgs) 


4.50 ± 0.50 


5.84 ± 0.08 


R(S.o) 


0.86 ± 0.07 


0.15 ±0.05 


reff(K) 


6000 ± 250 


/ 


^(Lg) 


0.87 ± 0.20 


/ 


Vo (mag) 


12.30 ±0.15 


/ 


My (mag) 


4.95 ±0.15 


/ 


d(pc) 


296 ± 12 


/ 



fitted by creating synthetic spectra in MOOG using the derived 
atmospheric parameters, which confirm the obtained metalicity. 
A detailed description of the method used t o determ i ne th e atmo- 
spheric parameters can be found in iPe Smedt et al.l (12012 . Sect. 
3). These atmospheric parameters correspond very well with the 
effective temperature and surface gravity derived from the SED 
fitting process. 

The atmospheric parameters derived from the spectral anal- 
ysis can be used to get a mass estimate for both the MS and 
sdB component independently from the SED fitting process, by 
using evolutionary models. Solar scaled evolutionary tracks of 
iBertelli et al.l (l2008l) with a metalicity of Z = 0.005 are compared 
to the derived surface gravity (logg - 4.5 + 0.5 dex) and effec- 
tive temperature (reif - 6000 + 250 K). The best fitting model 
has a mass of Mms - 0.86 + 0.15 Mq, with as corresponding age 
12 + 5 Gyr Using the mass ratio from the radial velocity curves, 
a mass of MsdB - 0.54 + 0.10 Mq is obtained for the sdB com- 
ponent. This corresponds within error to the canonical value of 
MsdB = 0.47 + 0.05 Mo. 

Using the newly obtained mass of the sdB component to de- 
rive its surface gravity from the gravitational redshift we find 
loggsdB = 5.84 + 0.08 dex. Furthermore, the radius, luminosity 
and magnitudes in the V band can be derived using the methods 
described in Sect.|5] the results are shown in Table[8] Most of the 
derived parameters are compatible with those derived in Sect.|5] 

7. Discussion and Conclusions 

From an analysis of literature photometry and observed spectra, 
detailed astrophysical parameters of PG 1 104+243 have been es- 
tablished. The surface gravity determined from the binary SED 
fitting method agrees very well with the surface gravity deter- 
mined from the observed gravitational redshift. Furthermore, the 
long time-base observations made it possible to accurately estab- 
lish the orbital period at 753 + 3 d. 

The sdB component is consistent with a canonical post-core- 
helium-flash model with a mass of -0.47 Mq, formed through 
stable Roche lobe overflow. This is supported by the essentially 
circular orbit as stable mass transfer is the only known process 
that can circularize a long-period binary system. 

When comparing our results with the distribution predicted 
by BPS models, we find that the orbital period of PG 1 104-H243 
is higher tha n the most likel y outcome at just over 100 days re- 
sulting from iHan et alJ (12003^ Fig. 21). However, the period is 
in agree ment with the observed estimate of three to four years 
noted by iGreen etaP (12001), and with BPS models computed 
within the y-formalism (Nelemans 2010). Our results therefore 
provide stron g observational cons traints for which BPS mod- 
els are viable. iClausen et an(l2012h published several limitations 



on the limiting mass ratio for dynamically stable mass trans- 
fer during the RLOF phase, the way that mass is lost to space 
and the transfer of orbital energy to the envelope during a com- 
mon envelope phase, based on the observed period distribution 
of sdB + MS binaries. They state that if long-period (Poit > 75 
d) sdB+MS binaries exist, mass transfer on the red giant branch 
(RGB) is stable, and if sdB + MS binaries with periods over 
250 d exists, mass that leaves the system carries away a spe- 
cific angular momentum proportional to that of the donor. In 
both cases the efficiency at which orbital energy is transferred to 
the envelope during the common envelope phase is below 75 %. 
Furthermore, in several of their BPS models the resulting dis- 
fi'ibutions for sdB + MS bina ries would include periods up to 
500 days ( IClausen et alj|20l"2l Fig. 16). However, this formation 
channel cannot explain the low mass of the cool companion. In 
this formation channel, the cool companion would be expected 
to have a mass ~ IMq. 

The low mass of the MS component (0.735 + 0.07 Mq) in- 
dicates that during the stable Roche-Lobe overflow phase only 
a small amount of the mass was accreted, and the major part 
must have been lost to infinity. The high evolutionary age of the 
MS component (~ 13000 Myr) is further evidence of this. If the 
companion had accreted a significant fraction of the original en- 
velope of the sdB component, its progenitor should have had a 
much lower mass as the current 0.735 + 0.07 Mq. In this case, 
the MS component should be practically un-evolved, as e.g. the 
0.6 Mq track of Fig. |6] which hardly evolves at all in a Hubble 
time. As the time since the mass transfer ended cannot be longer 
than the core-helium burning lifetime of only 100 Myr, no sig- 
nificant evolution of the MS component can have taken place 
sin ce then. 

iDeca et aP (l2012h published a study of the long-period 
sdB-i-K system PG 1018-047, which shares many similar prop- 
erties with PG 1104-1-243. It has a period of 760 + 6 d, a mass 
ratio of Msab/M^ - 0.63 + 0.11, and they find an effective 
temperature of 30500 + 200 K and logg of 5.50 + 0.02 dex 
for the sdB component. However, while the cool companion 
in PG 1104+243 contributes -52% of the Hght in the V-band, 
in PG 1018-047 it only contributes ~6%, and is photometri- 
cally and spectro scopically consistent with a mid-K type star. 
iDeca et al.l(l2012l) speculate that the orbit of PG 1018-047 may 
be quite eccentric, and that it can have f ormed through the 
merger scenario of IClausen & Wad3 (1201 1 *). meaning that the 
K-star was not involved in the evolution of the sdB star. The 
extremely low eccentricity of PG 1 104+243, makes such a sce- 
nario very unlikely, and we conclude that it is the first sdB+MS 
system that shows consistent evidence for being formed through 
post-Roche-lobe overflow. 

Thus, we have used high-resolution spectroscopy to solve 
the orbits of both components in an sdB+MS binary system, and 
have for the first time derived accurate and consistent physical 
parameters for both components. Furthermore, the accurate ra- 
dial velocities allowed the first measurement of gravitational red- 
shift in an sdB binary. PG 1 104+243 is part of an ongoing long- 
term observing program of sdB+MS binaries with HERMES at 
Mercator, and an analysis of the complete sample will make it 
possible to refine several essential parameters in the current for- 
mation channels. 
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